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The centrality dependence of the midrapidity charged-particle multiplicity density (|r/|<l) is pre- 
sented for Au+Au and Cu+Cu collisions at RHIC over a broad range of collision energies. The 
multiplicity measured in the Cu+Cu system is found to be similar to that measured in the Au+Au 
system, for an equivalent -/Vp art , with the observed factorization in energy and centrality still per- 
sistent in the smaller Cu+Cu system. The extent of the similarities observed for bulk particle 
production is tested by a comparative analysis of the inclusive transverse momentum distributions 
for Au+Au and Cu+Cu collisions near midrapidity. It is found that, within the uncertainties of 
the data, the ratio of yields between the various energies for both Au+Au and Cu+Cu systems are 
similar and constant with centrality, both in the bulk yields as well as a function of p T , up to at 
least 4 GeV/c. The effects of multiple nucleon collisions that strongly increase with centrality and 
energy appear to only play a minor role in bulk and intermediate transverse momentum particle 
production. 

PACS numbers: 25.75.-q,25.75.Dw 



The experimental program at the Relativistic Heavy- 
Ion Collider (RHIC) at Brookhaven National Laboratory 
has opened a new era in the study of QCD matter. Mea- 
surements of the midrapidity density of charged particles 
have, from the start of the RHIC program, provided an 
estimate of the energy density created in these collisions 
as well as an important baseline against which to test core 
features of particle production in models. The wealth of 
data on bulk charged-particle production collected in the 
first 5 years of RHIC operations has revealed the presence 
of various features in the data that appear to be uni- 
versal, including a factorization of the charged-particle 
production into two components, one that depends on 
collision energy, or ^/s NN , and the other on the central- 
ity of the colliding heavy- ions [l[. The energy depen- 
dence of the midrapidity charged-particle density appears 
to be consistent with a logarithmic rise with ^/s NN Q. 
This logarithmic dependence resembles that in elemen- 
tary p+p and e + + e _ collisions [3d and was not a priori 
expected for relativistic heavy-ion collisions. In partic- 
ular, initial expectations of a large increase in particle 



production with collision energy for central heavy-ion col- 
lisions due to multiple binary nucleon scatterings and the 
increased formation of mini-jets were not realized. Fur- 
ther studies revealed a relatively weak dependence of the 
produced midrapidity charged-particle multiplicity den- 
sity per participant pair on the collision centrality. From 
a two-component picture [J] , the part of the multiplicity 
which scales with the number of binary collisions (N co n) 
contributes at the level of only «13%, across all energies 
studied so far at RHIC 0,0. 

This paper reports results from RHIC for the midra- 
pidity charged-particle density in Cu+Cu collisions at 
y/s NN = 22.4, 62.4 and 200 GeV as a function of cen- 
trality. The data correspond to the top 50% of the total 
inelastic cross-section, and range from an average num- 
ber of participating nucleons (N part ) = 20 to 108. The 
previously published study of the centrality dependence 
of Au+Au multiplicity also covered the top 50% of the to- 
tal inelastic cross-section, corresponding to (iV part ) > 60. 
Thus, the full investigation of bulk charged-particle pro- 
duction as a function of iV part was limited by the wide gap 
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between the results of p+p collisions and the lowest mea- 
sured Au+Au centrality bin. For Au+Au, this left room 
for different scenarios of particle production evolution in 
heavy-ions, and made it impossible to distinguish various 
evolution scenarios (for example, a smooth rise or a sud- 
den jump at low A par t values). One way to circumvent 
this is via a detailed mapping of the particle production 
as a function of A part , from values corresponding to mid- 
central Au+Au collisions all the way down to those in 
d+Au collisions. The Cu+Cu data is therefore ideally 
suited for this study as the average number of partici- 
pants in the most central data class ((A part ) max = 108) 
overlaps with the peripheral Au+Au measurements, and 
the average number of participants in the most peripheral 
bin ( (Ap ar t ) mm = 20) reaches down to the upper limit for 
central d+Au collisions. In addition, for the same num- 
ber of participants as in Au+Au collisions, the relative 
uncertainty in the fractional cross-section of Cu+Cu col- 
lisions is smaller, enabling a more precise measurement. 
Also discussed in this paper, a subsequent re-analysis 
of the Au+Au data, using the techniques developed for 
Cu+Cu data analysis, extends the overlap in A part down 
to ~20 participants. 

The new Cu+Cu collision data were recorded by the 
PHOBOS experiment during the 2005 RHIC run. The 
primary hardware triggering was performed, as in prior 
PHOBOS heavy-ion data by two sets of 16 scintilla- 
tor "Paddle" counter arrays that cover a pseudorapidity 
of 3.2 < \tj\ < 4.5 (77 = -fn(tan(0/2))). The read-out 
from the full PHOBOS detector was initiated by the oc- 
currence of one or more scintillator hits on each array 
within At < 10 ns. This trigger prevents the loss of cen- 
tral collision events that may occur for Cu+Cu collisions 
if the primary hardware trigger is based solely on the 
zero-degree calorimeters (ZDC's), which detect spectator 
neutrons. For this analysis, the energy and time informa- 
tion from the ZDC's were used as an additional hardware 
cross-check to remove any contamination from beam-gas 
background events for mid-central collisions. The resul- 
tant data were then passed through a series of off-line 
quality checks before undertaking the physics analysis. 
The two most important of these checks include a cut 
on events with a time difference, between the Paddles, of 
5 ns or less to filter out beam-gas collisions, and a further 
cut on events that had a valid reconstructed primary ver- 
tex, determined as detailed below, within 10 cm of the 
nominal collision position. 

The primary algorithms for vertex reconstruction in 
PHOBOS were initially optimized for high multiplicity, 
central Au+Au collisions. For the Cu+Cu data, where 
multiplicity on average is much lower, the nominal re- 
construction procedures are found to be inefficient, even 
for mid-central data. In order to maximize the usability 
of this dataset, two new vertex reconstruction techniques 
were designed to allow for the analysis of the more periph- 
eral events. The first method finds the vertex position 



from hits in the single-layer Octagon detector (|7y| < 3.2) 
using a probabilistic approach applied to the energy de- 
posited at each hit Q. The overall vertex resolution at- 
tained by this technique is found to be a ~ 0.6 — 1.0 cm, 
dependent on the multiplicity of a given event. In order 
to obtain a better resolution along the beam axis (z), a 
complementary method based on all possible two-point 
straight tracks is used. Straight line trajectories are as- 
sumed as the magnetic field in the region of the vertex 
detector or the inner six planes of the Spectrometer (used 
in the vertex reconstruction) is negligible [8j . If the two- 
point track, formed from the Vertex detector and the first 
two silicon layers of the Spectrometer, passes close to the 
(x, y) orbit position of the beam, then the z-coordinate 
at the orbit position is considered a candidate for the 
collision vertex. From this ensemble of tracks, the z- 
position is determined from the most frequently found 
position. As a quality assurance measure, it is required 
that these two (independent) vertex positions are within 
|Az| < 3 cm of each other. This requirement removes 
all mis-identified collision vertices without detriment to 
the efficiency. Agreement to within 3 cm represents a 3<r 
acceptance of the lower resolution Octagon-based ver- 
tex. The overall vertex resolution for this procedure is 
determined by the two-point track vertex reconstruction 
method and is found to be a < 0.1 cm for high multi- 
plicity events and a ~ 0.25 cm for low multiplicity data. 
This overall improvement in vertex reconstruction effi- 
ciency, while still maintaining a good resolution, has al- 
lowed all track-based analyses (including those reported 
here) to be carried out in the Cu+Cu collision system. 

To provide a measure of the centrality of the colli- 
sion, and subsequent assignment of model-based param- 
eters such as the number of participants, the data are 
divided into bins of fractional cross-section. These bins 
allow for the comparison of data between collision sys- 
tems and energies. The first step toward determining the 
cross-section binning is to estimate the trigger and vertex 
reconstruction efficiency. Two methods were employed 
for the Cu+Cu data. The "Paddle hit" method counts 
the number of scintillator Paddle slats hit and compares 
them to the number from a fully simulated Monte Carlo 
(MC) sample of events. This method was also used in the 
analysis of 62.4-200 GeV Au+Au data. The inefficiency 
is estimated from the loss of events in the region of lowest 
number of Paddle slats hit, see for example Ref. [l[. A 
second technique, employed for both 19.6 GeV Au+Au 
and 200 GeV d+Au collisions, known as the shape match- 
ing method, compares the distribution of hit multiplic- 
ity in the Octagon detector to that in the MC simula- 
tion [||. The MC distribution is scaled to match the 
maximum multiplicity in data (assuming 100% efficiency 
at maximum multiplicity) such that a ratio of the mul- 
tiplicity distributions in data and in MC corresponds to 
the efficiency. Both of these techniques were used with 
two different MC generators (hijing [§] and AMPT (Iol |) 
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FIG. 1: The charged-particle pseudorapidity density at midrapidity (\rj\ < 1), scaled by the number of participating nucleon 
pairs, for (a) Au+Au collisions and (b) Cu+Cu collisions. The shaded ellipses in panel (a) represent the uncertainty on the 
new peripheral data, unshaded ellipses are from Ref. [l5|]. I n panel (b), the lines represent the Au+Au data (no errors) for 
comparison, where the 19.6 GeV data is scaled by 1.055 (derived from (151 ]') to account for the difference in collision energy of 
22.4 GeV for the Cu+Cu data. The open symbols at -/Vp art = 2 represent the inelastic p+p data at equivalent energies to the 
data, as in Ref. @]. The ellipses represent 90% C.L. systematic errors. 



and GEANT |ll| simulations of the detector response to 
determine the efficiency. The combined trigger and ver- 
tex reconstruction efficiencies for Cu+Cu collisions were 
found to be 84+5%, 75+5%, and 79+5% for 200, 62.4 
and 22.4 GeV respectively. A small difference (~5% in 
determined efficiency) was found between methods, and 
this was considered in the final systematic uncertainty of 
the centrality determination. For this analysis, the cen- 
trality cross-section bin selection was based on the total 
energy deposited in the Octagon detector. As a cross- 
check of auto-correlation biases, the centrality bin selec- 
tion was also determined using the sum of energy in the 
Paddles, truncated to the 24 (of 32) lowest energy signal 
slats to remove Landau-tail fluctuations. No difference in 
the results between the two centrality methods was found 
for peripheral to mid-central data. A 3% difference was 
found only for the most central data point, which is cov- 
ered by the systematic uncertainty on the measurement. 
The average number of participating nucleons, (-/V part ), 
for a given centrality bin was extracted from the initial 
Glauber-based values in the HIJING MC simulation for 
the associated centrality cuts on the final fully simulated 
Paddle and Octagon centrality distributions. System- 
atic uncertainties on (N palt } were determined using dif- 



ferent event generators including the uncertainties in the 
efficiency of both triggering and vertex reconstruction. 
For more details on the PHOBOS centrality methods see 
Refs. [H,[I3. 

The midrapidity charged-particle pseudorapidity den- 
sity measurement in this report is based on hit patterns 
in the dual-plane Vertex silicon detector ||. Two recon- 
structed hits in either the top or bottom pair of planes, 
form a "tracklet". Among all possible hit combinations, 
or tracklets, in these two detector planes, only those that 
point back to the pre-determined vertex position (within 
a distance of closest approach of 6 mm) are utilized for 
the final measurement, greatly reducing the number of 
fake tracklet candidates. The charged-particle multiplic- 
ity is determined by counting the number of tracklet can- 
didates and applying a correction factor to account for 
the geometrical acceptance of the detector, reconstruc- 
tion efficiency, and vertex resolution [||. Possible effects 
on the tracklet data, due to the changing vertex resolu- 
tion with multiplicity, were investigated with no system- 
atic effects observed. The results were further checked for 
possible auto-correlation biases (due to the new vertex 
reconstruction procedures) via comparison to the p T in- 
tegrated invariant yields [13| and the single-layer analysis 
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FIG. 2: The ratio between two collision energy datasets (200 GeV / X GeV) of the charged particle pseudorapidity densities 
from Fig. [1] Panel (a) shows the data ratios from Au+Au collisions 15]. Panel (b), shows the Cu+Cu data, along with lines 
representing the Au+Au data (no errors) for comparison. In panel (b) the 19.6 GeV Au+Au data is scaled by 1.055 to account 
for the difference in collision energy between the low energy Au+Au (19.6 GeV) and Cu+Cu (22.4 GeV) data. The triangles 
(at iVp ar t = 2) represent the inelastic p+p data at equivalent energies to the Au+Au data, as in Ref. Q. Vertical error bars 
are combined statistical and systematic 1-a uncertainties. 



using the Octagon detector [l4j. All results are found to 
be consistent. 

The charged-particle pseudorapidity density, at midra- 
pidity, scaled by the number of participating nucleon 
pairs is shown in Fig. Q] for Au+Au and Cu+Cu colli- 
sions. Panel (a) shows the published Au+Au data [l5j 
and for each energy four new peripheral data points, with 
shaded error ellipses. The new Au+Au data (analyzed 
using the techniques developed for Cu+Cu) extend the 
studies of charged particles at mid-rapidity from a frac- 
tional cross-section of 50% to a peripheral value of 70%. 
The triggering and vertex reconstruction efficiency are 
100% in this region. With these new peripheral data it is 
now apparent that the charged-particle multiplicities at 
midrapidity tend towards the inelastic p+p data. In ad- 
dition, the new peripheral Au+Au data provide a large 
overlap in (N pait ) with the Cu+Cu data, Fig. QJi. The 
Au+Au data points are also shown in Fig. [TJd as a line, 
in comparison to the Cu+Cu data. The charged particle 
yields are higher in Cu+Cu collisions than for inelastic 
p+p data, and, as for Au+Au collisions, the multiplicity 
per participant pair increases with centrality and collision 
energy. Within the systematic uncertainty, the trends in 
the centrality dependence of the charged particle produc- 



tion are the same across systems. The degree to which 
the trends are the same is somewhat surprising in the 
context that the number of participants differ by a factor 
of three for the same cross-sectional fraction, due to the 
much smaller size of the Cu nucleus. 

A more precise way to look at trends in the data is 
to take ratios of the midrapidity charged particle yields 
at different collision energies, as a function of centrality. 
In this way a large fraction of the systematic uncertain- 
ties in the measurement cancel, as detailed in leav- 
ing a combined (statistical and systematic) uncertainty 
at the level of a few percent. The primary reason for 
this cancellation is that the analysis uses the same de- 
tector, reconstruction technique, triggering requirements 
and centrality determination method across the different 
data samples. Systematic uncertainties arising from a 
common source are thus neutralized. The result of this 
analysis for all PHOBOS data on Au+Au and Cu+Cu 
collisions is shown in Fig. [2] These ratios are taken as 
a function of fractional cross-section and shown versus 
the corresponding average iV part value between the two 
energies. The small differences in (iV part ) for the same 
fraction of cross-section for different collision energies has 
a negligible effect within a given collision system. Per- 
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forming the analysis in the Au+Au system for precisely 
matched (N pait ) values yielded the same results [6|. Fur- 
thermore, the ratios obtained using different centrality 
methods are consistent within the shown uncertainties. 

We observe two striking features in the data. First, 
even within the reduced uncertainties, the slope of the ra- 
tios are consistent with zero for both Cu+Cu and Au+Au 
systems as a function of centrality. This observation in 
the data clearly illustrates that the factorization in cen- 
trality and collision energy of the midrapidity charged- 
particle yields is a feature of both the large Au+Au sys- 
tem as well as the smaller Cu+Cu system. As the ra- 
tios show no measurable centrality dependence, due to 
the factorization into energy and centrality components, 
the trends of dN^/dr]]^^ with centrality seen in Fig. Q] 
are not driven by the collision energy and thus must be 
a consequence of the collision geometry, i.e. the volume 
of the interaction as defined by the number of partici- 
pating nucleons. Second, we observe that the increase 
in charged-particle yields with collision energy is similar 
between the Cu+Cu and Au+Au systems. In particu- 
lar, the ratio of midrapidity charged-particle yields from 
200 and 62.4 GeV are the same in Cu+Cu and Au+Au 
systems, within the systematic uncertainties. The offset 
in the ratio of 200/19.6 Au+Au compared to 200/22.4 
Cu+Cu is largely understandable in terms of the dif- 
ference in the lower collision energies (i.e. 19.6 versus 
22.4 GeV). 

The extent of the similarities in the charged-particle 
ratios for Cu+Cu and Au+Au collisions can be fur- 
ther tested by a comparative analysis of the inclu- 
sive charged-particle transverse momentum (p T ) distri- 
butions. For this, data from Ref. [13] are used to form a 
ratio of charged-particle yields between y/s NN = 200 and 
62.4 GeV as function of centrality and p T . The results, 
shown in Fig. [3] for three bins of p T , are independent of 
centrality, in agreement with those from the bulk yields, 
even for transverse momenta up to p T ^4 GeV/c. One 
may expect that the produced yields at intermediate- to 
high-p T would be predominantly formed from hard par- 
tonic collisions, i.e. those more likely to produce jets. As 
the hard collisions should follow a scaling with N co \\ then 
a distinct N pait dependence of the ratio should be ob- 
served for different p T ranges. In this regard, it is sur- 
prising to find no evidence of a centrality dependence 
apparent in the data at high-p T . The increase of the av- 
erage ratio value with p T is expected due to the harder 
spectra at 200 GeV compared to 62.4 GeV, as also ob- 
served in p+p collisions. A more definitive statement 
at high-p T would require data with higher statistics and 
lower systematic uncertainties. 

It is important to point out that the effects observed 
in the p T dependence favor a volume-dominated scenario 
of particle production, and it will be interesting to see 
whether the further differential measurements (particle 
species dependent effects) support this conclusion. Com- 
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FIG. 3: Ratio (200 GeV / 62.4 GeV) of charged-particle yields 
measured at three different (p T ) values versus centrality in 
Au+Au (closed symbols) and Cu+Cu (open) collisions. Ver- 
tical error bars represent the statistical uncertainty in the 
measurement, brackets represent 90% C.L. systematic errors. 
The data are from Refs. QJ, [TJ, \vj§. 



parisons of proton to pion ratios in Au+Au collisions at 
200 and 62.4 GeV are consistent with our observations 
for unidentified charged-particles [l8j]- From our study, 
it is clear that effects of multiple-nucleon collisions which 
strongly grow with centrality and energy play only a mi- 
nor role, if any, in bulk charged-particle production. 

In conclusion, this paper reports new data for the 
midrapidity charged-particle density in Cu+Cu and pe- 
ripheral Au+Au collisions at RHIC. The charged parti- 
cle density per participant pair is found to be similar be- 
tween peripheral Au+Au and Cu+Cu, for the same num- 
ber of participants. The midrapidity yields of the smaller 
Cu+Cu system appear to factorize into collision energy 
and centrality dependent components, in much the same 
manner as seen in the Au+Au system. The observed 
centrality independence of the ratios, between different 
center-of-mass energies, further illustrate that the colli- 
sion geometry drives the bulk yield of charged particle 
production at RHIC. This apparent general feature also 
appears to hold as a function of transverse momentum 
for p T up to at least ~ 4 GeV/c. 
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